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Abstract 

We study the PT-symmetric boundary conditions for "spin" -related 6- 
interactions and the corresponding integrabihty for both bosonic and 
fermionic many-body systems. The spectra and bound states are dis- 
cussed in detail for spin-| particle systems. 

Key words: Spin-coupling, 5-interaction, PT-symmetry, Integrability 
PACS number(s): 02.30.Ik, ll.30.Er, 03.65.Fd 

Integrable models play significant roles in statistical and condensed matter physics. Many of 
these models can be exactly solved in terms of an algebraic or coordinate "Bethe Ansatz method" 
[1], see e.g., [2, 3] for spin chain and ladder models and [4, 5] for (continuous variable) quantum 
mechanical models. In particular, the quantum mechanical solvable models describing a particle 
moving in a local singular potential concentrated at one or a discrete number of points have been 
extensively discussed, see e.g. [6, 7, 8] and references therein. One dimensional problems with 
point interactions at, say, the origin (x = 0) can be characterized by the boundary conditions 
imposed on the wave function at x = 0, which is equivalent to two particles with contact 
interactions. The integrability of one dimensional quantum mechanical many-body problems 
with general contact interactions between two particles, has been studied in [9] and the bound 
states and scattering matrices are calculated for both bosonic and fermionic statistics. The 
results are generalized to the case of quantum mechanical systems with "spin" -related contact 
interactions, namely, the boundary conditions describing the contact interactions are dependent 
on the spin states of the particles [10]. 

Recently, the complex generalization of conventional quantum mechanics has been investi- 
gated [11]. In stead of the standard formulation of quantum mechanics in terms of Hermitian 
Hamiltonians, quantum mechanical models with space-time reflection symmetry (PT symmetry) 
have been constructed and studied for continuous interaction potentials [12]. For point interac- 
tion potentials a systematic description of the boundary conditions and the spectra properties 
for self-adjoint, PT-symmetric systems and systems with real spectra have been presented, and 
the corresponding integrability of one dimensional many body systems with these kinds of point 



interactions are studied [13]. In this article we study the PT-symmetric boundary conditions 
for "spin" -related (5-interactions and the corresponding integrability for many-body case. 



The family of the usual point interactions for the one dimensional Schrodingcr operator 
— -^ can be described by 2 x 2 matrices. In particular, the boundary conditions describing the 
(^-potential interactions have the following form 

n \ / ,o \ 

(1) 





0+ \ / \ / 0- 
where c G R, (p{x) is the scalar wave function of two spinless particles with relative coordinate 
X. (1) also describes two particles with spin s but without any spin coupling between the 
particles when they meet (i.e. x = 0), in this case tp represents any one of the components of 
the wave function. It is easily verified that the boundary condition (1) is both self-adjoint and 
PT-symmetric. 

Now consider two particles with spin s and have both J-interactions and spin couplings when 
they meet. For a particle with spin s, the wave function has n = 2s + 1 components. Therefore 
two particles with 5-interactions have a general boundary condition, described in the center of 
mass coordinate system: 

where the wave function ij; and its derivative if)' are n^-dimensional column vectors, C is an 
X V? matrix, and I2 is the v? x identity matrix. 

For general C the interactions described by the boundary condition (2) is neither self- 
adjoint nor PT-symmetric. From the symmetry condition of the Schrodinger operator, < 
>L2{R,c") - < ^'-^^ >L2(R,c")= 0) for any n, d G C°°(M \ {0}), one has that 
(2) is self-adjoint when C = C\ where f stands for the conjugate and transpose. 

Applying PT operator to (2) we have 
Hence the boundary conditions for the function PTil) are given by 






which coincides with (2) if and only if 





(3) 



This implies that C* = C, i.e., C is a real matrix. 
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Therefore the boundary condition (2) is self-adjoint when C is Hermitian and PT-symmetric 
when C is real. When C is real symmetric, then the interaction described by the boundary 
condition (2) is both self-adjoint and PT-symmetric. When C = CI2, (2) is reduced to the 
case (1) without spin-coupling interactions when two particles meet. In the following we study 
quantum systems with contact interactions described by the boundary condition (2). We first 
consider two spin-s particles case. The Hamiltonian is then of the form 

^ = ^"5"^^^^ + ^^^^'^^"'^^^' 

where h is an x matrix. If the matrix h is proportional to the unit matrix I2, H is reduced 
to the usual two-particle Hamiltonian with contact interactions but no spin coupling. 

Let Ca, a = 1, ...,n, be the basis (column) vector with the a-th component as 1 and the rest 
components 0. The wave function of the system (4) is of the form 

n 

^ = 4'ap{xi,X2)ea^ep. (5) 

In the center of mass coordinate system, X = {xi + X2)/2, x = Xi — X2, the operator (4) has the 
form 

The functions (j) = from the domain of this operator satisfy the following boundary 

conditions at x = 0, 

<?!'a/3(0+, X) = (l)ap{0~,X), a,P = 1, n, 

where the indices of the matrix h are arranged as 11, 12, In; 21, 22, 2n; nl, n2, ...,nn. 
(7) is just the boundary condition (2) with C = h. Here h acts on the basis vectors of the two 

n 

particles as h ep) = ^ hafi,-y\&^ ® ex. 

7,A=1 

According to the statistics ip is symmetric (resp. antisymmetric) under the interchange of 
the two particles if s is an integer (resp. half integer). Let k\ and k2 be the momenta of the 
two particles. In the region xi < X2, in terms of Bethe hypothesis the wave function has the 
following form 

where U12 and U21 are x 1 column matrices. In the region xi > X2, 

where according to the symmetry or antisymmetry conditions, P^^ = p^"^ for bosons and P^"^ = 
—p^^ for fermions, being the operator on the ra^ x 1 column that interchanges the spins of 
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the two particles. Substituting (8) and (9) into the boundary conditions (7), we get 

\ U12 + U21 = P^'^{ui2+U2l), 

\ iki2{u21 - U12) = hP^^{ui2+U2l) +iki2P'^^{ui2 - U2l), 

where ki2 = {ki — k2)/2. EUminating the term P^^u\2 from (10) we obtain the relation 

U2l = Yi^Ul2 , (11) 

where 

Yi^ = [2ifci2l2 - h]-\2iki2P^^ + h]. (12) 

For a system of N identical particles with (5-interactions, the Hamiltonian is given by 

N q2 n 

H = - JZ^'^N + Y^hijS{xi-Xj), (13) 
i=i ^-^i i<j 

where Iat is the x identity matrix, hij is an operator acting on the i-th and j-th bases 
as h and the rest as identity, e.g., hi2 = ^ (8) I3 <8) ...Iat, with Ij the n x n identity matrix acting 
on the i-th. particle. The wave function in a given region, say xi < X2 < ■■■ < xn, is of the form 

n 

* = X] 4>ai,...,aN{xi,...,XN)eai<Si ...<Siea^ 

ai,...,ajv=l , . 

= U12 jve*^'^^^!"^'^^^^''"""''"'^^^^^ + U21 jYe*^''^^!"'"'^^^^"'""""'"'^^^^^ ^ ' 

+{N\ — 2) other terms, 



where kj, j = 1, N, are the momentum of the j-th particle, u are n x 1 matrices. The wave 
functions in the other regions are determined from (14) by the requirement of symmetry (for 
bosons) or antisymmetry (for fermions). Along any plane Xi = Xj+i, i G 1, 2, A'' — 1, we have 

^aiQ2.--ajQ:j+i...ajv ~ ^^^.^a^ ^^aia2...aj+iQ:j...ajv ' (f^) 

where 

Y^:^, = ['^ika,a,+jN - hjj+i]-\2ik^.^.^,P^^+^ + hjj+i]. (16) 

Here kajaj+i = {kaj — kaj+i)/^ play the role of momenta and pj^+^ = p>^'^^ for bosons and 
pjj+i — _pjj+i fQj. fermions, where p'^'^^ is the operator on the x 1 column u that inter- 
changes the spins of particles j and j + 1. 

For consistency the operators Y in (16) must satisfy the Yang-Baxter equation with spectral 
parameter [4, 5], 

Y-rn,rn+l-irrn+l,rn+2-y-rn,m+l -i^m+l,m+2-i^m,m+l-i^m+l,m+2 

ij kj ki ki kj ij ' 



or 



ij kj ki ki kj ij 



h 



(19) 



if m, r, s are all unequal, and 

ij ji ^ kl hi "ij ^ ' 

if m,r,s,q are all unequal. By a straightforward calculation it can be shown that the operator 
Y given by (16) satisfies all the Yang-Baxter relations if 

[hij,P'^]=0. (18) 

Therefore if the Hamiltonian operators for the spin coupling commute with the spin permutation 
operator, the A^-body quantum system (13) can be exactly solved. The wave function is then 
given by (14) and (15) with the energy E = J2iLi kf. 

For the case of spin-^, ^ is a 4 x 4 matrix. A matrix satisfying (18) is generally of the form 

( a e\ e\ c \ 

e2 f g 63 

62 5 / 63 

yd 64 64 b J 

where a,b,c, d, f,g, 61,62, €3,64 G C. When 62 = e^, 64 = 63, d = c* and a,b,f,g G M, the 
interaction is reduced to the self-adjoint case [10]. when a, b, c, d, /, g, ei, 62, 63, 64 G R, the 
interaction is PT-symmetric, and when 62 = ei, 64 = 63, d = c and a,b,c, f,g,6i,63 G M, the 
interaction is both self-adjoint and PT-symmetric. 

One can compare this kind of quantum mechanical integrable systems with the integrable spin 
chain models. Let V denote a complex vector space. A matrix R taking values in EnddV V) 
is called a solution of the quantum Yang-Baxter equation (QYBE) without spectral parameters, 
if it satisfies 7^-12^13^23 = T^2'iT^ViT^i2^ where IZij denotes the matrix on the complex vector 
space V ® V ® V , acting as R on the i-th and the j-th components and as identity on the 
other components. When V is two dimensional, the solutions of QYBE include the ones of the 
form (19) [14], although generally (19) is not a solution of the QYBE. A solution R in principle 
corresponds to a kind of spin coupling operator between the nearest neighbor spins and gives 
rise to an integrable spin chain model by Baxterization [15]. Therefore for an A^-body system 
to be integrable, the spin coupling in the contact interaction (19) is allowed to be more general 
than the spin coupling in a Heisenberg spin chain model with nearest neighbors interactions. 

The spectra of (13) with h given by (19) are generally complex. They are real when the 
spectra of h are real. The self-adjoint case is included in the parameter family that gives real 
spectra. The real spectra related to the PT-symmetric case are divided into two parts: self- 
adjoint and non-Hermitian. 

When h gives purely real spectra, The many-body system (13) possesses also bound states. 
For N = 2, from (10) the bound states have the form. 



C+Q.A 



V'^ = n„e^l^^-^^1, a = l,...,n^ (20) 
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where Ua is the common a-th eigenvector of the operators h and P , with eigenvalue Aq, 
such that hua = A^Uq. and c + aA^ < 0, P^^Ua = Ua- The eigenvalue of the Hamitonian H 
corresponding to the bound state (20) is — (c + akaf /2. In stead of a unique bound state for 
(5-interaction with boundary condition (1), here we have ■n? bound states. For the A?^-particle 
system the bound states are then given by 

tpa=Vae 2 ^r>j' a = l,...,n , (21) 

where Va is the wave function of the spin part. It can be checked that satisfy the boundary 
condition (7) at Xi = xj for any i / j G l,...,N. The spin wave function v here satisfies 
P^-'va = Va and hijVa = ^aVa, for any i / j. The energy of the bound state tp^ is -Eq = 
— (c + aAa)^N{N^ — 1)/12. The corresponding scattering matrices can be similarly studied for 
real spectra case. 

We have studied the boundary conditions with spin-coupling (^-interactions including self- 
adjoint, PT-symmetric and real spectrum families. The corresponding integrability for both 
bosonic and fermionic many-body systems are investigated. Detailed examples are given to 
the case of spin-^ particle systems, where the spectra and bound states are also discussed. 
The results can be straightforwardly generalized to higher spin particle systems. The contact 
interaction we considered here is just (5-interaction. For general contact interactions with spin- 
couplings, the PT-symmetries, spectra problems, and integrability of many-body problem could 
be much more complicated and remain to be investigated further. 
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